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Abstract 

We review the theories of a few quantum phase transitions in two-dimensional corre- 
lated electron systems and discuss their application to the cuprate high temperature 
superconductors. The coupled-ladder antiferromagnet displays a transition between 
the Neel state and a spin gap paramagnet with a sharp 5 = 1 exciton: we develop a 
careful argument which eventually establishes that this transition is described by the 
familiar 0(3) tp 4 field theory in 2+1 dimensions. Crucial to this argument is the role 
played by the quantum Berry phases. We illustrate this role in a one-dimensional 
example where our results can tested against those known from bosonization. On 
bipartite lattices in two dimensions, we study the influence of Berry phases on the 
quantum transition involving the loss of Neel order, and show that they induce 
bond-centered charge order (e.g. spin Peierls order) in the paramagnetic phase. We 
extend this theory of magnetic transitions in insulators to that between an ordinary 
d-w&ve superconductor and one with co-existing spin-density-wave order. Finally, 
we discuss quantum transitions between superconductors involving changes in the 
Cooper pair wavefunction, as in a transition between d x 2_ y 2 and d x 2_ y 2 + id xy su- 
perconductors. A mean-field theory for this transition is provided by the usual BCS 
theory; however, BCS theory fails near the critical point, and we present the re- 
quired field-theoretic extension. Our discussion includes a perspective on the phase 
diagram of the cuprate superconductors, in the context of recent experiments and 
the transitions discussed here. 
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1 Introduction 



The theory of quantum phase transitions has emerged in the last decade as 
a powerful tool for the description of the unconventional electronic proper- 
ties of a variety of correlated systems. This article will provide an elementary 
overview of some theoretical models of quantum phase transitions in two spa- 
tial dimensions, along with a brief discussion of their experimental applica- 
tions. One reason for our focus on two dimensions is, of course, our interest 
in the properties of the cuprate superconductors, whose interesting electronic 
excitations reside within a Cu02 square lattice. However, not coincidentally, 
two is also the spatial dimension in which many of the most natural class of 
non-trivial quantum critical points reside; three is the "upper-critical dimen- 
sion" of many critical points, at and above which the physical properties are 
adequately described by a mean field theory (possibly extended with some low 
order fluctuation corrections). 

This article is addressed to newcomers to the study of quantum phase transi- 
tions who desire an introduction to some of the experimentally relevant theo- 
retical questions under current discussion. We will assume the reader is well- 
versed in the standard theory of classical critical phenomena, and so our focus 
will be on novel aspects of quantum phase transitions, including the role of 
quantum Berry phases and of excitations with fermionic statistics. Earlier re- 
views by this author [1,2] provide detailed discussions of background material 
which we will occasionally refer the reader to; most of our discussion here goes 
beyond these reviews, and focuses on current theoretical and experimental de- 
velopments. 

We begin in Section 2 by describing the magnetic quantum transition in the 
coupled ladder antiferromagnet. High precision numerical studies of the criti- 
cal properties of this transition are now available, and these serve as a useful 
check on the theoretical models. We map the critical fluctuations of the mag- 
netic Neel order on to the 0(3) tp A field theory in 2+1 dimensions, along with 
additional Berry phase terms which reflect the commutation relations of the 
individual quantum spins. We illustrate the role of the Berry phases by first 
considering a one-dimensional antiferromagnet in Section 3, where the results 
can be compared with those obtained by bosonization methods. We return to 
our discussion of two-dimensional quantum spin models in Section 4, where 
we extend the methods developed in Section 3: we demonstrate that the Berry 
phases induce bond-centered charge order in the paramagnetic phase of an- 
tiferomagnets with full square lattice symmetry, while they can be safely ne- 
glected in the coupled ladder antiferromagnet. The remaining sections extend 
these theories to doped antiferromagnets with mobile charge carriers. Section 5 
considers magnetic transitions in superconductors, while Section 6 describes 
transitions between BCS superconductors with differing pair wavefunctions. 
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Fig. 1. The coupled ladder antiferromagnet. Spins (S = 1/2) are placed on the 
sites, the A links are shown as full lines, and the B links as dashed lines. 

2 Coupled ladder antiferromagnet 



We will begin our discussion by describing the quantum phase transition in 
a simple two-dimensional model of antiferromagnetically coupled S = 1/2 
Heisenberg spins. At the microscopic level, this model does not describe the 
spin fluctuations in the cuprate superconductors. However, a fairly subtle ar- 
gument, which we develop in the following sections and culminate in Section 5, 
shows that the universal properties of the simple critical point described in this 
section are identical to those of a magnetic ordering transition in the cuprate 
superconductors. The model presented here also provides a useful description 
of other transition metal oxides with a spin gap [3]. 

We consider the "ladder" Hamiltonian 

H e = J ^ • s, + aj Yl s* • Si (1) 



where Sj are spin- 1/2 operators on the sites of the coupled-ladder lattice shown 
in Fig 1, with the A links forming decoupled two- leg ladders while the B links 
couple the ladders as shown. The ground state of H depends only on the 
dimensionless coupling A, and we will describe the low temperature properties 
as a function of A. We will restrict our attention to J > and < A < 1. 

We will begin with a physical discussion of the phases and excitations of the 
coupled ladder antiferromagnet, Hi in Section 2.1. We will propose a quantum 
field-theoretical description of this model in Section 2.2: we will verify that 
the limiting regimes of the field theory contain excitations whose quantum 
numbers are in accord with the phases discussed in Section 2.1, and will then 
use the field theory to describe the critical point. 
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Fig. 2. Schematic of the quantum paramagnet ground state for small A. The ovals 
represent singlet valence bond pairs. 

2.1 Phases and their excitations 

Let us first consider the case where A is close to 1. Exactly at A = 1, H is 
identical to the square lattice Heisenberg antiferromagnet, and this is known to 
have long-range, magnetic Neel order in its ground state i.e. the spin-rotation 
symmetry is broken and the spins have a non-zero, staggered, expectation 
value in the ground state with 

(S,) = Vj N n, (2) 



where n is some fixed unit vector in spin space, and rjj is ±1 on the two 
sublattices. This long-range order is expected to be preserved for a finite range 
of A close to 1. The low-lying excitations above the ground state consist of slow 
spatial deformations in the orientation n: these are the familiar spin waves, 
and they can carry arbitrarily low energy i.e. the phase is 'gapless'. There 
are two polarizations of spin waves at each wavevector k = (k x , k y ) (measured 
from the antiferromagnetic ordering wavevector), and they have excitation 
energy e k = {(? x k\ + c^) 1//2 , with c x ,c y the spin-wave velocities in the two 
spatial directions. 

Let us turn now to the vicinity of A = 0. Exactly at A = 0, H is the Hamilto- 
nian of a set of decoupled spin ladders. Such spin ladders are known to have 
a paramagnetic ground state, with spin rotation symmetry preserved, and an 
energy gap to all excitations [4] . A caricature of the ground state is sketched 
in Fig 2: spins on opposite rungs of the ladder pair in valence bond singlets in 
a manner which preserves all lattice symmetries. Excitations are now formed 
by breaking a valence bond, which leads to a three-fold degenerate state with 
total spin S = 1, as shown in Fig 3; this broken bond can hop from site-to-site, 
leading to a triplet quasiparticle excitation. Note that this quasiparticle is not 
a spin-wave (or equivalently, a 'magnon') but is more properly referred to as 
a spin 1 exciton. For A small, but not exactly 0, we expect that the ground 
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Fig. 3. (a) Cartoon picture of the bosonic 5 = 1 exciton of the paramagnet. (b) 
Fission of the S = 1 excitation into two 5=1/2 spinons. The spinons are connected 
by a "string" of valence bonds (denoted by dashed ovals) which are not able to 
resonate with their environment, and some of which lie on weaker bonds; this string 
costs a finite energy per unit length and leads to the confinement of spinons. 

state will remain a gapped paramagnet, and the spin exciton will now move 
in two dimensions. We parameterize its energy at small wavevectors by 



Su = A H „ ■ 

2A 



(3) 



where A is the spin gap, and c x , c y are velocities. Fig 3 also presents a simple 
argument which shows that the S = 1 exciton cannot fission into two S =1/2 
'spinons'. 

The very distinct symmetry signatures of the ground states and excitations 
between A ~ 1 and A ~ make it clear that the two limits cannot be contin- 
uously connected. It is known [5,6] that there is an intermediate second-order 
phase transition at A = A c ~ 0.3. Both the spin gap A and the Neel order 
parameter Nq vanish continuously as A c is approached from either side. 



2.2 Quantum field theory 



A convenient field-theoretic description of the two phases introduced in Sec- 
tion 2.1, and of the critical point which separates them, is obtained by the 
coherent state path integral. A detailed derivation of this approach has been 
provided by the author in Chapter 13 of Ref. [1] which the reader should now 
consult; we will only highlight a few important points here. 

The coherent state path integral is best introduced by considering first a single 
spin S with angular momentum S and Hamiltonian 

H = -h ■ S (4) 
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where h is a static magnetic field. The (2S+1) eigenvalues of Hq are, of course, 
— m|h|, where m = —S, —S+ 1, . . . , S — 1, S, and so we also know its partition 
function Z at a temperature T. The coherent state path integral represents 
Z in terms of an integral over all possible closed curves, N(r), on the surface 
of a unit sphere, each curve representing a possible history of the precessing 
spin in imaginary time, r: 



Z = j VN(t)5(N 2 (t) - 1) exp (-iS j A T (N(r))dT + J drSh • N(r) 



.(5) 



The first term is the so-called 'Berry phase', and A t (N(t))cIt is defined to be 
the oriented area of the spherical triangle defined by N(r), NfY + rfr), and an 
arbitrary reference point No (which is usually chosen to be the north pole). 
We will see that A T behaves in many respects like the time-component of a 
gauge field, and this accounts for the suggestive notation. All physical results 
should be independent of the choice of No, and it is easy to see that changes 
in N amount to gauge transformations of A T . So if we change N to N , then 
the resulting A' T is related to A T by 

A' T = A T - <9 r 0(r) (6) 



where 0(r) measures the oriented area of the spherical triangle defined by 
N(r), N , and N . 

We now apply the above path integral independently to every site of the 
coupled ladder antiferromagnet: we obtain a path integral over fields N,(t), 
where j is a site index as in (1). As discussed in Chapter 13 of Ref. [1], it 
is convenient to now take a spatial continuum limit, to allow focus on the 
important low energy excitations. This proceeds by the parameterization 

a 2 



N,-(r) = i&nfa.rtyl " (a 2 /S) 2 L 2 (r„T) + -L(r„r), (7) 



where rjj = ±1 on the two subattices, a is the lattice spacing, r = (x,y), 
n(r, r) is a unit length continuum field representing the staggered Neel order 
parameter, and L(r, r) is the continuum field for the conserved, uniform, mag- 
netization density (the normalization chosen in (7) is that required to make 
/ d 2 rL(r, t) exactly equal to the total spin). Next, a gradient expansion in 
performed on the effective action for Hi, and after integrating out the L fluc- 
tuations we obtain the following path integral for the partition function of the 
coupled ladder antiferromagnet: 



Z e = J Vn(r,r)5(n 2 (r,T) - 1) exp -iSj^Vj J drAr{n( rj , r)) 
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I d 2 rdr ((d T n) 2 + c 2 x (d x n) 2 + c 2 y (d y n) 2 ) , (8) 



where the velocities c x = JSay(l + A) (3 + A), c y = JSa^J 2(3 + A), and the 

coupling constant g = (2a/S)[(3 + A) 2 /(2 + 2A)] 1//4 . We will now view g as the 
tuning parameter, and values of g on either side of a critical value g c place 
the system in one of the two phases discussed in Section 2.1; notice that g is 
a decreasing function of A for < A < 1, and so the "strong" coupling phase 
g > g c corresponds to the spin gap phase with A < A c , while g < g c is the 
A > A c Neel phase. 

The first term in (8) is a residual Berry phase, and this remains expressed 
on the underlying lattice: there is no natural way to take its continuum limit. 
Notice that this term induces complex weights in the partition function and so 
it has no analog in the theory of classical phase transitions. This Berry phase 
plays a central role in determining the phases and critical points of many 
quantum spin systems, and describing its consequences is one of the central 
purposes of this article. However, our discussions in the following sections will 
show that at least for the case of the coupled ladder antiferromagnet the Berry 
phase in (8) can be safely neglected, and so we will drop it in the remainder 
of this subsection. 

Upon omission of the Berry phase in (8), Zi becomes a standard path integral, 
familiar in well established theories of classical phase transitions in 3 = 2 + 1 
dimensions. For instance, after interpreting r as a third dimension (i.e. r — > z), 
and n as the local magnetization in a three-dimensional ferromagnet (e.g. 
iron), then is simply its classical partition function at a finite 'temperature' 
T c i, with the coupling constant g ~ T c i. From the theory of such ferromagnets 
it is known that a convenient description of the physical properties of Zt 
emerges not in a formulation in terms of the fixed length field n (the 'hard 
spin' representation), but in a formulation using a 'soft spin' field </? Q (r, r) 
(a = x,y,z) with no length constraint. We connect the two formulations by 
a coarse-graining transformation, so that tp a represents the average of n over 
some small coarse-graining volume of spacetime. The new form of Zi is now 



where keep in mind that the Berry phase in (8) has been dropped, s is now 
the new tuning parameter (replacing g in (8)), and the quartic interaction u 
has replaced the fixed length constraint in (8). The two phases of Section 2.1 
are now separated by a critical value s = s c , with the A < A c spin gap phase 



Z e = 



J V<p a (r, r) exp - J d 2 rdr {(d rVa ) 2 + c 2 x (d x ip a ) 2 + c 2 y (d yVa ) 2 
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corresponding to s > s c , while the A > A c Neel phase appears for s < s c . 



The well-known properties of the model (9), when translated into the quantum 
framework of interest here, are easily seen to reproduce the physical properties 
of the phases discussed in Section 2.1. 

For s < s c , Zi in (9) has a "ferromagnetic" phase associated with the breaking 
of spin rotation invariance, and the development of the expectation value 
(ip a ) = 5 az N ] to lowest order in u, the spontaneous magnetization is N = 



Qs/u and is determined by the minimum of the "Mexican hat" potential 
for the tp a field. After using (7) this "ferromagnetic" order in the classical 
model is seen to correspond to the antiferromagnetic Neel order (2) in the 
quantum model. A standard Gaussian fluctuation analysis about this ordered 
states shows that susceptiblity transverse to the spontaneous moment is given 
by 

X±(k,u n ) = (\<f x (k,uj n )\ 2 ) = (\ Vy (k,u n )\ 2 ) = 1 (10) 

w n c x^x ' ^V^y 



where uj n is a frequency associated with imaginary time. After analytically 
continuing to real frequencies, u, this yields the quantum response function 
to an applied staggered field transverse to the Neel order: 

^•»>-&>+U->*- (11> 



This response function has a pole, representing the two spin-wave excitations, 

1/2 

at energy (c 2 J) . This is in complete accord with our discussion of 

the properties of the Neel phase in Section 2.1. 

For s > s c , the "disordered" phase of Z e in (9) corresponds to the quantum 
paramagnet with a spin gap. Here, a computation of the susceptibility to 
lowest order in u yields 

* ( *-"> = A' +C re + C ,>j- M " (12) 



where we have already analytically continued to real frequencies, A = y/s, 
and (for now) the residue Z — 1. This response function is independent of 
the direction of the applied staggered field, and so its pole represents the 

5 = 1 exciton at energy (A 2 + (? x k 2 x + Cyk^j . This is again in accord with 
the discussion of the spin gap phase in Section 2.1: the exciton dispersion 
(3) is obtained from the present result in a small momentum expansion. The 
pole structure in (12) actually holds to all orders in an expansion in u [1]: at 
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quantum temperature T = 0, such corrections renormalize the values of the 
spin gap A and the quasiparticle residue 2, and induce additional spectral 
weight associated with multiple excitonic continua at energies larger than 
3A, but they do not broaden the quasiparticle pole. Physically, this happens 
because the exciton is the lowest energy excitation with spin S = 1, and so 
conservation of total spin prevents implies that this excitation has an infinite 
lifetime. Only at T > does a finite lifetime appear, associated with scattering 
off the thermally excited density of pre-exiting excitons. 

The present field theoretical framework can also be extended to obtain a the- 
ory for the critical point at s = s c (A = A c ), as has been discussed in some 
detail in [1]. In its application to the classical ferromagnet, the critical prop- 
erties of (9) have been understood in some detail in studies of the ip A field 
theory, and accurate values of the critical exponents have been obtained. In a 
recent quantum Monte Carlo study of the coupled ladder antiferromagnet Hg, 
Matsumoto et al. [7] determined the critical exponents at A = A c to impres- 
sive accuracy, and all values are in accord with the exponents of the <p A field 
theory. This agreement is independent numerical support for the omission of 
the Berry phases in proceeding from (8) to (9); separate theoretical arguments 
for the same approximation will emerge in the following sections. As we ap- 
proach the critical point with s \ s c , the pole structure in (12) continues to 
hold at low frequencies, but the spin gap vanishes as A ~ (s — s c ) zu and the 
residue vanishes as Z ~ (s — s c ) vl/ , where z — 1, v, and rj are standard critical 
exponents. Right at s = s c there are no quasiparticle excitions, and the dy- 
namic susceptiblity instead displays a branch-cut representing the continuum 
of critical excitations: 



For more information on the physics of this critical continuum, refer to [1]. 



3 Quantum XY chain 

This section will divert from the main subject of this paper by considering some 
simple examples of quantum phase transitions in one dimension. Our purpose 
here is to illustrate the consequences of the Berry phases in (8) by considering 
a simplified model for which they can be evaluated completely; we will also 
be able to test our method by comparing the results with those obtained 
by a bosonization analysis of the same one-dimensional model. We will then 
proceed to apply closely related methods in two dimensions in Section 4. 
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We will evaluate the generalization of (8) in one dimension for case in which 
the Hamiltonian has on a global U(l) spin symmetry, and the anisotropies are 
such that the spins prefer to lie within the x-y plane in spin space. In this case 
we can parameterize n by a single angle 9: 

n = (cos 0, sin0, 0); (14) 



with these simplifications, we will find that powerful duality methods enable 
us to obtain considerable insight into the physics a partition function with 
Berry phases and complex weights. With the parameterization (14) and the 
restriction to one spatial dimension, we expect the resulting Zg to provide a 
description of a large class of spin S quantum antiferromagnets in one dimen- 
sion; to be specfic, we expect it will model the Hamiltonian 

Hxy — J\ ^ (S X jS x j + i + S yj Syj +1 + (S zj S z j +1 ) 
j 

+ J2^2 (SxjSx,j+2 + SyjSy,j+2 + (S z jS z j +2 ) , (15) 
3 

where > are first/second neighbor antiferromagnetic exchange constants 
(J 2 is not too large), and the anisotropy \(\ < 1 induces a preference for the 
spins to reside in the x-y plane. We can reasonably expect that the phase dia- 
gram of Hxy as a function of ( and J2 / J\ is similar to that of Zg as a function 
of g. The former is already known by an earlier bosonization analysis [9] and 
will provide an instructive comparison with the latter. 

Our computation begins by a transformation of the Berry phase in (8) to 
an alternative form in one dimension. We discretize spacetime into a square 
lattice of sites, j, (j is now a spacetime index, unlike the spatial index in (15), 
and its interpretation should be clear from the context). For this spacetime 
lattice we define Ajp, with \i = x,r, to be the oriented area of the spherical 
triangle formed by rij, rij +/i , and a fixed reference N . Then, as illustrated in 
Fig 4, we can use the identity 

X) VjAjr = fy^A^Aju (16) 

3 3 

to relate the Berry phase to the flux piercing the plaquettes of the space- 
time lattice; here £j is a fixed integer field on the sites, J, of the dual lattice, 
and a convenient choice is £j — (1 — (— 1)^)/2 i.e. £j is unity on every odd 
column, and zero on every even column. The symbol A M represents a discrete 
lattice derivative, with A^fj = fj + /x — fj. 

Now we use the in-plane form of n in (14) to our advantage. The flux e^A^Aju 
through any plaquette is clearly equal to the surface area of the quadrilateral 
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Fig. 4. Diagrammatic proof of the identity (16) which relates the Berry phase to 
the sum over the flux in every plaquette in the odd columns. 

on the unit sphere formed by the \\j on the four vertices of the plaquette. For 
planar spins, such a quadrilateral lies entirely on the equator, and its area 
is generically zero. The only exception occurs when a vortex in 9 pierces the 
plaquette: in this case, the quadrilateral is the entire equator and the flux is 
the area of a hemisphere = ±27r. So we have established that the Berry phases 
in (8) generate the factor 

exp —i2irS l ? x (vortex number of plaquette j) (17) 



for easy-plane spin systems in one dimension. As the vortex number is always 
an integer, we see that the factor in (17) is always unity for integer S, and the 
Berry phases can be ignored for these values of S. 

The remaining terms in the action of (8) can also be simplified for the present 
situation. On the discrete spacetime lattice, and with (14), these terms become 
—(l/g)YljCOs(A fl 9j), which is the usual action of a "classical" XY model in 
2 = 1 + 1 dimensions. As is standard, duality manipulations on this XY 
model are more conveniently performed in a periodic Gaussian, 'Villain', rep- 
resentation of the weights of the model. In this manner, we finally obtain the 
generalization of the partition function (8) to U(l) symmetric, quantum spin 
S chains: 

Zxy = S / II de i ex P ( ~o~ XXA^j - 2?rm iM ) 2 - i2nS ^ W m i 

Wd 3 V 19 31* D 



where m JAt are integers on the links of the direct lattice which, when summed, 
produce the periodic Gaussian weights for A^Oj, the □ represents a sum over 
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plaquettes, and e^A^m^ are the vortex numbers of the plaquettes. The ex- 
pression (18) is amenable to an exact duality transformation which produces 
a representation in which all weights are positive: we apply the Poisson sum- 
mation identity to the sum over on each link, and replace it by a sum 
over integers Jj M : 



Here, and henceforth, we drop overall normalization constants in front of par- 
tition functions. Now the integral over the 9j can be performed independently 
on each site, and only imposes the constraint A M Jj M = 0. We solve this con- 
straint by writing Jj M = e^ u A u pj, with the pj integers on the sites of the dual 
lattice. So we obtain the exact mapping of (18) to [8] 



where the "heights" Hj on the sites of the dual lattice are restricted to the 
values 



The expression (20) is a canonical height (or 'solid-on-solid') model in two 
dimensions. The only influence of the quantum Berry phase are the restrictions 
on the allowed values of the heights of the "interface" : for integer S the heights 
take all integer values just as in the standard solid-on-solid model, while for 
half-integer S, the heights are restricted to half-odd-integer (integer) values 
on the odd (even) columns of the square lattice. As we have made a number 
of approximations in the deriving Zxy from the path integral of spin systems, 
we shall be interested in a generalized class of height models with additional 
short-range non-local couplings, apart from the simple terms appearing in 
(20). With these additional couplings, standard methods from the theory of 
interface models can be used to show [8] that the generalized height models 
describe phases not only in the easy-plane, small \(\, limit of (15), but also 
those in the easy- axis, large \(\ limit when the spins prefer to lie in the ±z 
directions. In all of this analysis it is essential that we not relax the restrictions 
on the allowed values of the heights, as these are fundamental consequences 
of the quantized value of the spin S. We describe the main results for two sets 
of S values in turn. 





(20) 



Hj ^ Pj ~ Si j 



(21) 
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3.1 S integer 



The standard model with integer heights on every site has two phases. For 
small g, the interface is rough; this is dual to the classical or quantum XY 
phase with quasi-long-range spin order in the x-y plane, with power-law decay 
of spin correlations. For the quantum spin chain, we may refer to this phase 
as a Tomonaga-Luttinger (TL) liquid, as it shares many of the characteristics 
of other TL liquids. 

For large g, the interface is smooth; this is dual to the "disordered" phase of 
the classical XY model. For the quantum spin chain, this large g phase has 
exponential decay of spin correlations in imaginary time, and therefore a spin 
gap. Its characteristics are rather similar to the spin gap phase discussed in 
Section 2: the primary excitation is a gapped S = 1 exciton with an infinite 
lifetime at low energies. As the spin anisotropy is reduced by increasing ( in 
(15), this spin gap phase connects continuously to the well-known Haldane 
gap state of integer spin chains with SU(2) symmetry. 

The transition between the small and large g phases is second-order, and 
is described by the standard Kosterlitz-Thouless theory of the classical XY 
model in two dimensions. 

3.2 S half- odd-integer 

The staggering of the heights between integer and half-odd-integer values in- 
troduces crucial differences in the phase diagram, and several new phases are 
allowed in the generalized phase diagram of such height models. However, it 
is believed [8] that the specific height model (20), with only a local gradient 
squared coupling, has only one rough phase for all values of g, and this corre- 
sponds, as above, to a TL phase of the quantum spin chain. The other phases 
appear upon allowing additional short-range couplings in (20), and these can 
reasonably be expected to arise from terms in the original Hamiltonian omit- 
ted in the duality mappings; as we noted above, in this generalization it is 
essential that we maintain the restrictions on the allowed height values. The 
new phases are associated with smooth interfaces, and as in Section 3.1 such 
phases have a spin gap. However, the staggering of the heights induces crucial 
differences in the structure of such spin gap states, as we describe below. 

A crucial characteristic of a smooth interface is that there is a well-defined 
value of the average height (Hj) in every pure state. Furthermore, for the 
staggered height models under consideration here, every value of the average 
(Hj) breaks a lattice translational symmetry of the quantum spin chain. For 
instance if (Hj) = (modulo integers) then the even columns of the square 
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Fig. 5. (a) Schematic of the state with bond-centered charge (BC) or spin-Peierls 
order for the S = 1/2 quantum spin chain. The ellipses represent valence bonds as 
in Fig 2, and these have spontaneously chosen one of two possible dimerizations of 
the chain, (b) S = 1/2 spin excitations of the BC state. Unlike the case in Fig 3b, 
these spinons are free to move far apart from each other: this is because we are in 
one dimension and the dimerization is spontaneous. Such spinons are confined in 
two dimensions whether the BC order is spontaneous or explicitly present in the 
Hamiltonian. 

lattice (which have integer allowed heights) have been picked out, and we 
may expect the average energy, (Sj x Sj + x >x + Sj y Sj + i ;y ) on the links associated 
with the even columns to be smaller (say) than those on the odd columns. 
Conversely, an interface with (Hj) = 1/2 (modulo integers) reverses the role 
of the odd and even columns. The choice between these two sets of average 
heights corresponds to a spontaneously broken Z 2 symmetry in the ground 
state, which now has 'spin-Peierls' or 'bond-centered charge' (BC) order. The 
nature of these ordered states, and of their excitations, is illustrated in Fig 5: 
they have a spin gap, but the lowest-lying excitations are fractionalized S = 
1/2 spinons which reside on the domain wall between the two BC states. 

Other possible values of (Hj) in the smooth interface phase lead to other 
discrete broken symmetries in the quantum spin chain. As we describe in 
Ref. [8], the values (Hj) = 1/4,3/4 (modulo integers) correspond to states 
with Ising spin order i.e. (Sj Z ) = ±r]jN . All other values of (H 3 ) correspond 
to four-fold degenerate states with co-existing Ising and BC order. 

We sketch the phase diagram of the model (20) for S half-odd-integer in Fig 6, 
as obtained by a bosonization analysis [9]. The TL, BC, and Ising phases are 
all separated by second-order transition lines, and the three lines meet at a 
single multi-critical point: all the critical properties can be described by an 
extension of the Kosterlitz-Thouless renormalization group flows. Analysis of 
the generalized height models discussed here leads to a phase diagram with 
an identical topology and critical properties [8]. 
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Fig. 6. Phase diagram of the 5 = 1/2 quantum spin chain Hxy as obtained by 
abelian bosonization [9]. Identical results are obtained from a generalized square 
lattice height model in which the heights take half-odd-integer (integer) values on 
the even (odd) columns. 

4 Berry phases in two dimensions 

We are now prepared to describe the evaluation of the 2+1 dimensional path 
integral with Berry phases in (8). As in the 1+1 dimensional case in Section 3 
the evaluation is performed on a discrete lattice in spacetime: the degrees of 
freedom are three component unit vectors, n,, on the sites, j of a cubic lattice. 
On such a lattice we can rewrite (8) as 

Z = [ II dn^n) - 1) ex P ( -7T £ n i • n i+A - iS Yl VjAjr ) , (22) 

J j V 9 3 J 

where the sum over \x extends over the three directions x,y,T. In the application 
to the coupled latter antiferromagnet, the spatial anisotropy and doubled unit 
cell of Hi led to an anisotropy in the spatial derivative terms in (8). In the 
discrete form (22) this reduced symmetry requires a corresponding variation in 
the couplings g associated with the couplings between the various sites on the 
cubic spacetime lattice: we have dropped this variation in (22) and so, strictly 
speaking, this partition function applies only to quantum antiferromagnets 
which have the full symmetry of the square lattice and have a single spin S on 
each site. An example of such an antiferromagnet is the so-called J1-J2 model 

H=J 1 J £S i -S j + J 2 £ S i -S i (23) 

(ij) nnn ij 

with exchange interaction J\ > between square lattice nearest neighbors 
(ij), and J 2 > between next-nearest-neighbors (nnn). As in Section 3, we 
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exclude large values of J2 and expect that the phases of H as a function of 
Jij J\ will be described by those of Z in (22) as a function of g. Later we will 
also discuss the consequences of doubling the unit cell and introducing spatial 
anisotropy in (22), which is required to obtain a description of the coupled 
ladder model in (1). 

For small g, we expect Z to be in the ordered Neel state, with broken spin 
rotation symmetry. It is not difficult to account for the Berry phases in a spin- 
wave fluctuation analysis about this ordered state, and they only have minor 
consequences [8]. 

The Berry phases are more crucial in the paramagnet state, with full spin 
rotation symmetry, expected at large g. Because this is a strongly coupled 
system in the 3 dimensions with a non-abelian symmetry, no exact duality 
analysis is possible (unlike the simpler situation in Section 3). We proceed 
with the strategy discussed in Ref. [8]: we replace the integral over the n,- 
by an integral over the Aj^ the variables that appear in the Berry phases. 
Formally, this can be done by introducing new 'dummy' variables A, M and 
rewriting (22) by introducing factors of unity on each link; this leads to 



Z = J n^exp l-iS^njAjA J HdnjSirf-lWAj^-A^) 

m \ 3 / 3 




In the first expression, if the integral over the A, M is performed first, we trivially 
return to (22); however in the second expression we perform the integral over 
the rij variables first at the cost of introducing an unknown effective action Sa 
for the Ajp. In principle, evaluation of Sa may be performed order-by-order in 
a "high temperature" expansion in 1/g: we match correlators of the Aj^ flux 
with those of the A 3IX flux evaluated in the integral over the n, with positive 
weights determined only by the 1/g term in (22). Rather than undertaking this 
laborious calculation, we can guess essential features of the effective action Sa 
from some general constraints. First, correlations in the rij decay exponentially 
rapidly for large g (with a correlation length ~ l/ln(g)), and so Sa should 
be local. Second, it should be invariant under the lattice form of the gauge 
transformation (6) 

A' m = A jlJt - A^j (25) 
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associated with the change in the reference point on the unit sphere from No 
to N , with <j)j equal to the area of the spherical triangle formed by n^, N 
and N . Finally the area of any triangle on the sphere is uncertain modulo 4-7T 
and so the effective action should be invariant under 

A m - A m + 4tt. (26) 



The simplest local action which is invariant under (25) and (26) is that of 
compact U(l) quantum electrodynamics and so we have 

Z = J II dA m ex P E cos (\^xKA jX ^j -iSj2 VjAj^j , (27) 



for large g; comparison with the large g expansion shows that the coupling 
e 2 ~ g 2 . As in Section 3 our analysis is aided by replacing the cosine interaction 
in (27) by a Villain sum over periodic Gaussians: 

E (l^uxKAjx - 27vq Jfl ) -iSj2 VjAjr I , 



z=E/n^exp(-^i;(i 



where q^ are integers on the links of the dual lattice, which pierce the pla- 
quettes of the direct lattice. 

We will now perform a series of exact manipulations on (28): our final result, in 
(35), will be a demonstration of its exact equivalence to another height model, 
but now with the interface 3=2+1 dimensional. As in Section 3, this interface 
model is obtained after a duality mapping and has only positive weights. This 
last fact, of course, makes it much more amenable to a standard statistical 
analysis. This first step in the duality transformation is to rewrite (28) by the 
Poisson summation formula: 



E ex P -^Efo^MjA - 2vrg 



2\ 



= E ex P ( E a \ - 1 E l^vxa^AvAjx ) , (29) 
Km) V 1 J D 1 J 

where (like q^) is an integer-valued vector field on the links of the dual 
lattice. Next, we write the Berry phase in a form more amenable to duality 
transformations. Choose a 'background' = a° flux which satisfies 

e^xKa°jx = Vj^r, (30) 



17 



t • • i • 
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x +1 x * x +1 * X 



Fig. 7. Specification of the non-zero values of the fixed field a^. The circles are the 
sites of the direct lattice, j, while the crosses are the sites of the dual lattice, J; the 
latter are also offset by half a lattice spacing in the direction out of the paper (the 
fj, = T direction). The are all zero for fj, = t,x, while the only non-zero values of 
djy are shown above. Notice that the a flux obeys (30). 



where j is the direct lattice site in the center of the plaquette defined by the 
curl on the left-hand-side. Any integer- valued solution of (30) is an acceptable 
choice for a2 , and a convenient choice is shown in Fig 7. Using (30) to rewrite 
the Berry phase in (28), applying (29), and shifting by the integer 2Sa2 ! , 
we obtain a new exact representation of Z in (28): 

Z = J2 [ UdAj^exp ( E(%* - 2Sa Jfl ) 2 - iJ2 \^u\a^A v A jX ] .(31) 



The integral over the A, M can be performed independently on each link, and 
its only consequence is the imposition of the constraint e^ u \A u aj\ = 0. We 
solve this constraint by writing as the gradient of a integer-valued height 
hj on the sites of the dual lattice, and so obtain 

Z = E ex P f E( A ,% " 2£aJJ 2 ] • (32) 
{hj} V 5./* / 



This is the 2+1 dimensional height model in almost its final form. Before we 
proceed to an analysis of (32), we pause for a paragraph to make some paren- 
thetic remarks on the application of the above approach to one dimension. 

Precisely the same approach as above can be applied to SU(2) symmetric 
spin chains which we can model by the obvious generalization of (22) to 1+1 
dimensions. In the large g limit, we proceed with the same set of duality 
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XXX 

+1/8 -1/8 +1/8 -1/8 
(b) 



Fig. 8. Specification of the non-zero values of the fixed fields (a) Xj and (b) 
introduced in (34). The notational conventions are as in Fig 7. Only the ji = r 
components of are non-zero, and these are shown in (b). 

mappings, and in place of (32) we obtain the partition function 

Zi = £<Wex P L e l^{a,-2Si^f ] . (33) 



The zero gradient constraint on aj (to be contrasted with the zero curl con- 
straint on in 2+1 dimensions) means that it takes a single fixed value 
through the entire system. For integer S, the value aj = S minimizes the 
action (its action per site is e 2 S 2 /2, to be contrasted with the two-fold de- 
generate BC states with aj = S — 1, S + 1 with action per site e 2 (S 2 + l)/2), 
and yields a 'featureless' state with no broken symmetry: this is, of course 
the Haldane gap state expected for SU(2) symmetric, integer S, spin chains, 
as discussed in Section 3.1. For half-odd-integer S, there is BC order, as the 
action in minimized by the two-fold degenerate values = S + 1/2, S — 1/2. 
All of these results in 1+1 dimensions are in complete accord with the results 
obtained in Section 3 for large g, and with the known properties of quantum 
spin chains. 

Now we return to the 2+1 dimensional case of interest in this section, and 
proceed with our interrupted analysis of (32). We can convert the "frustra- 
tion" cij^ in (32) into offsets for the allowed height values (as in Section 3) by 
decomposing it into curl and divergence free parts and writing it in terms of 
new fixed fields, Xj and as follows: 

a Jll = A li X s +e IJa ,xA v y j x. (34) 



The values of these new fields are shown in Fig 8. Inserting (34) into (32), we 
can now write the height model in its simplest form [10-12] 



(35) 
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where 



Hj = hj - 2SXj (36) 

is the new height variable we shall work with. Notice that the have dropped 
out, while the Xj act only as fractional offsets (for S not an even integer) to 
the integer heights. From (36) we see that for half-odd-integer S the height is 
restricted to be an integer on one of the four sublattices, an integer plus 1/4 
on the second, an integer plus 1/2 on the third, and an integer plus 3/4 on the 
fourth; the fractional parts of these heights are as shown in Fig 8a; the steps 
between neighboring heights are always an integer plus 1/4, or an integer plus 
3/4. For S an odd integer, the heights are integers on one square sublattice, 
and half-odd-integers on the second sublattice. Finally for even integer S the 
offset has no effect and the height is an integer on all sites. We discuss these 
classes of S values in turn in the following subsections. 

(It is worth emphasizing here that the height model in (35) only offers a 
description of the large g paramagnetic phase of the model Z in (22) in 2+1 
dimensions. The small g Neel ordered state is not one of the phases of Z h , nor 
is the critical point describing the onset of magnetic order; an extension of Zh 
to include the Neel state was discussed in [8], but will not be presented here. 
This should be contrasted with the situation in 1+1 dimensions discussed in 
Section 3; there height models like Z X y in (20) can describe all the phases 
and critical points of the quantum spin chain, including the TL phase, which 
is a state with quasi-long-range Neel order). 



4-1 S even integer 

Unlike the two-dimensional case, three-dimensional height models generically 
have no roughening transition, and the interface is always smooth. With all 
heights integers, the smooth phase breaks no symmetries. So like the case of 
integer S in quantum spin chains (Section 3.1), square lattice antiferromagnets 
with S even integer can have a paramagnetic ground state with a spin gap and 
no broken symmetries. This is in accord with the exact ground state for a S = 2 
antiferromagnet on the square lattice found by Affleck et al. [14], the AKLT 
state. The structure of this state is identical to that of the paramagnetic state 
of the ip A field theory in 2+1 dimensions (Eqn (9) with c x = c y ), with a stable 
S = 1 spin exciton excitation. So the critical point between the magnetic Neel 
state and the spin gap state should also be described by the y? 4 field theory. 
As we noted above, this critical point is not contained within Zh- 
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Fig. 9. Mapping between the quantum dimer model and the interface model Z in 
(35). Each dimer on the direct lattice is associated with a step in height of ±3/4 
on the link of the dual lattice that crosses it. All other height steps are ±1/4. Each 
dimer represents a singlet valence between the sites, as in Fig 2. 

4-2 S half- odd-integer 



Now the heights of the interface model can take four possible values, which 
are integers plus the offsets on the four square sublattices shown in Fig 8a. 
As in Section 4.1, the interface is always smooth, but as in Section 3.2, any 
smooth interface must break a lattice symmetry with the development of bond- 
centered charge (BC) order and this allows a number of distinct spin gap 
ground states of the lattice antiferro magnet. Consistent with these theoretical 
predictions, numerical studies of the frustrated spin model (23) have produced 
evidence for BC order in the state contiguous to the Neel state: see [15,16] for 
recent results. 



It is useful, first, to obtain a simple physical interpretation of the interface 
model in the language of the S — 1/2 antiferromagnet [13]. From Fig 8a it 
is clear that nearest neighbor heights can differ either by 1/4 or 3/4 (modulo 
integers). To minimize the action in (35), we should choose the interface with 
the largest possible number of steps of ±1/4. However, the interface is frus- 
trated, and it is not possible to make all steps ±1/4 and at least a quarter 
of the steps must be ±3/4. Indeed, there is a precise one-to-one mapping be- 
tween interfaces with the minimal number of ±3/4 steps (we regard interfaces 
differing by a uniform integer shift in all heights as equivalent) and the dimer 
coverings of the square lattice: the proof of this claim is illustrated in Fig 9. We 
identify each dimer with a singlet valence bond between the spins (the ellipses 
in Fig 2), and so each interface corresponds to a quantum state with each 
spin locked in the a singlet valence bond with a particular nearest neighbor. 
Fluctuations of the interface in imaginary time between such configurations 
correspond to quantum tunneling events between such dimer states, and an 
effective Hamiltonian for this is provided by the quantum dimer model [17,18]. 

The nature of the possible smooth phases of the interface model are easy 
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(a) (b) 

Fig. 10. Sketch of the two simplest possible states with bond-centered charge (BC) 
order for S = 1/2 on the square lattice: (a) the columnar spin-Peierls states, and (b) 
plaquette state. The different values of the (Sj • S,-) on the links are encoded by the 
different line styles. Both states are 4-fold degenerate; an 8-fold degenerate state, 
with superposition of the above orders, also appears as a possible ground state of 
the generalized interface model. 

to determine from the above picture, and by standard techniques from sta- 
tistical theory which have been reviewed elsewhere [11,8,13,19]. Interfaces 
with average height (Hj) = 1/8,3/8,5/8,7/8 (modulo integers) correspond 
to the four-fold degenerate BC ordered states in Fig 10a, while those with 
(Hj) = 0, 1/4, 1/2, 3/4 (modulo integers) correspond to the four-fold degener- 
ate plaquette BC states in Fig 10b. All other values of (Hj) are associated with 
eight-fold degenerate BC states with a superposition of the orders in Fig 10a 
and b. 

It is interesting to note here that the symmetry of the state in Fig 10a is 
identical to the symmetry of the coupled ladder Hamiltonian illustrated in 
Fig 1: so the system has spontaneously decided to arrange itself in the form 
of a coupled ladder. Whether the BC ordering is spontaneous (as it is here) 
or explicitly present in the Hamiltonian (as in Section 2), the nature of the 
non-zero spin excitations should be the same. So the confinement argument of 
Fig 3b applies here too, and there are no free S = 1/2 spinons in the BC state. 
The elementary excitation is a S — 1 spin exciton. This should be contrasted 
with the S = 1/2 quantum spin chain in Section 3.2, where the BC state had 
spinon excitations. 

We are now finally in a position to explain a central step in the analysis of 
Section 2: the neglect of Berry phases in the theory of the magnetic critical 
point for the coupled ladder antiferromagnet. We have seen here that the 
primary effect of the Berry phases is to induce BC order in the spin gap 
state. However, such "order" is already present explicitly in the Hamiltonian 
of the coupled ladder antiferromagnet, which has two sites per unit cell: the 
Hamiltonian picks out one of the four BC states and "dimer" fluctuations to 
the other BC states are suppressed. So we need only focus on the fluctuations 
of the magnetic order, and these are captured by the (p 4 field theory. 
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The nature of the crossover between the Neel and BC states for S — 1/2 
quantum spin models with the full square lattice symmetry (like H in (23)) 
is far more subtle, and was recently discussed at some length in Ref [8]. Here, 
both the Neel and BC orders can become critical, and so we cannot naively 
assert the irrelevance of the Berry phases. One appealing scenario [15,8], sup- 
ported by recent numerical studies [15], is that the crossover occurs via two 
transitions: with increasing g, there is first a transition at g — g c \ to a state 
with co-existing Neel and BC order (only the BC fluctuations are critical at 
g c i are these are described by a Z 4 clock model), and then a transition at 
9 — 9c2 > 9ci where Neel order disappears. The second transition at g c2 is just 
as in the coupled ladder antiferromagnet, and is therefore described by the (p A 
field theory: there is well-formed BC order on the either side of g — g&, and 
it is immaterial here whether this order was explicitly present in the Hamilto- 
nian or spontaneously induced. Indeed, this was secretly the reason we chose 
to consider the coupled ladder antiferromagnet in our early discussion. 

4-3 S odd integer 

This case is similar to that S half-odd-integer, and we will not consider it in 
detail. The Berry phases again induce BC order in the spin gap state, but this 
order need only lead to a two-fold degeneracy. 



5 Magnetic transitions in d-wacve superconductors 

We have so far considered a variety of magnetic transitions in correlated 
"Mott" insulators, which have an electron density commensurate with the 
number of available orbitals. Cuprate high temperature superconductivity ap- 
pears when such a Mott insulator on the square lattice is doped by a finite 
density of mobile carriers. All of the available experimental evidence supports 
the proposal that, apart from some inhomogeneous charge ordered insulating 
states, the ground state of the doped Mott insulator is a superconductor. In 
particular, there is no clear-cut evidence for a metallic Fermi liquid state up 
to a moderate density of holes. Consequently, we will focus our attention on 
the superconducting ground state in the remainder of this paper. This is to 
be contrasted with the approach of [20] which examines transitions in metallic 
states at non-zero temperature. 

It is now natural to study the interplay of the Neel magnetic order found in 
the insulator and the superconducting order of the doped system. (The charge 
order discussed in Section 4 also plays an important role in the superconduc- 
tor, but we will not focus attention on this aspect here — see [21] for a recent 
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review.) The experimental evidence indicates that this crossover proceeds via 
an intermediate state with co- existence of magnetic and superconducting or- 
der. Furthermore, the magnetic order is no longer in the simple two-sublattice 
Neel configuration, but oscillates at a wavevector K ^ (n, n). So instead of 
expressing the electron spin as Sj a = r]jip a (rj), which is implied by (7) and 
the subsequent coarse-graining of n to ip a , we now write 



where & a (r,r) a 3-component, complex quantum field; $ a is required to be 
complex as along as the magnetic state does not have a two sublattice struc- 
ture, which is only the case for ordering wavevectors (ir, ir) and (0,7r). Also, 
the experimental evidence also indicates that the magnetic order is collinear, 
and not spiral; this is the case if e aj a 7 ($ /3 )( ( l ) *) = 0. 

With this reasoning, we are led to consider the quantum transition between 
a superconducting state with co-existing magnetic spin-density-wave order at 
wavevector K 7^ in the SC+SDW state) to an ordinary superconduc- 

tor = in the SC state). Clearly, the field § a should be a degree of 

freedom in any theory for this transition. However, in a superconductor there 
are additional fermionic excitations which carry spin S = 1/2, the Bogoliubov 
quasiparticles, and we have to examine their influence on the critical prop- 
erties. This is especially the case in a rf-wave superconductor, in which the 
energy of the quasiparticles vanishes at four 'nodal' points in the Brillouin 
zone. 

To address this issue, let us review the standard BCS mean-field theory for a 
d-w&ve superconductor on the square lattice. We consider the Hamiltonian 



where Cj a is the annihilation operator for an electron on site j with spin 
o~ =|, I, Cka is its Fourier transform to momentum space, Sk is the dispersion 
of the electrons (it is conventional to choose Sk = —2ti(cos(k x ) + cos(k y )) — 
2t 2 (cos(A; x -|- k y ) + cos(k x — k y )) — fi, with t 12 the first/second neighbor hopping 
and \i the chemical potential), and the J\ term is the same as that in (23) 
with 



and a a the Pauli matrices. We will consider the consequences of the J2 term in 
(23) for the superconductor in Section 6 below. Applying the BCS mean-field 




(37) 




(38) 




(39) 
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decoupling to H t j we obtain the Bogoliubov Hamiltonian 

Hbcs = E ^A a c kc - y E A m ( C it c i+A,l - c )i c )+m) + h - c - ( 40 ) 

k j/i 



For a wide range of parameters, the ground state energy optimized by a d x 2_ y 2 
wavefunction for the Cooper pairs: this corresponds to the choice A x = —A y = 
A x 2_ y 2. The value of A x 2_ y 2 is determined by minimizing the energy of the 
BCS state 

r d 2 k 

E BCS = Ji\A x 2_ y 2\ 2 -J—[E k - e k ] (41) 



where the fermionic quasiparticle dispersion is 
Eh 



' 2 21 1 / 2 

e k + | JiA a; 2_j / 2(cos k x — cos k y )\ . (42) 



The energy of the quasiparticles, Ek, vanishes at the four points (±Q,±Q) 
at which e k = 0. We are especially interested in the low energy quasiparticles 
in the vicinity of these points, and so we perform a gradient expansion of 
Hbcs near each of them. We label the points Q\ = (Q,Q), Qi = {—Q,Q), 
Qs = (-Q, -Q), Qa = (Q, -Q) and write 

c ja = fia(xj)e^ + h*{x 3 Y Q ^ + h„i ■>■;,)< + / 4CT (^)e^, (43) 



while assuming the fi-4, c (x) are slowly varying functions of x. We also intro- 
duce the bispinors *i = (/i T , f^, f n , -/ 3 f T ), and ^ 2 = (fih / 4 \, /21, -fl\), and 
then express Hbcs in terms of ^1,2 while performing a spatial gradient expan- 
sion. This yields the following effective action for the fermionic quasiparticles: 

= J drd 2 x (d T - iv F T z d x - iv A r x d y ) *i 

+ (d T ~ lV F T Z 8 y - IV A T X 8 X ) * 2 ] , (44) 

where we have rotated our spatial co-ordinate system by 45°, and the t x,z are 
4x4 matrices which are block diagonal, the blocks consisting of 2 x 2 Pauli 
matrices. The velocities vp^ A are given by the conical structure of E k near the 
Qi- 4 : we have v F = \V k e k \ k=Qa \ and v A = | JiA a .2_ J/ 2v / 2sin(Q)|. 

We are now in a position to discuss the quantum field theory of the transition 
from the SC+SDW to the SC state in term of the <& a and ^2 degrees of 
freedom. First, we can write down an action for $ a along the lines of that 
for ip a in the insulator: this will have a structure very similar to that of (9) 
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(the complex nature of $ Q induces additional quartic terms which we will 
not discuss here — the reader is referred to another review by the author [21] 
for details on this point). The action will also contain the terms in in 
(44), and it remains to discuss terms which may couple the $ Q and ^2. The 
simplest possible terms are cubic interaction terms like <3> Q \[>i\I/i etc. However 
by comparing the momentum dependencies in (37) and (43) with (39) it is 
easy to see that such terms will survive after averaging over space only if 
the special commensurability conditions K = (2Q,2Q), or K = (2(5,0), are 
satisfied. In general this will not be the case, and experimental evidence also 
does not support this possibility. Therefore, we will assume in the remainder 
of this section that these commensurability conditions are not satisfied: the 
theory for the case K = (2Q,2Q), Q = it /2 was considered in Ref. [22], and 
has some formal similarities to the theory to be considered in Section 6 in a 
different physical context. 

In the absence of special commensurability conditions, the simplest term which 
can couple the <3> a and ^1,2 is 



However a simple scaling argument shows that the coupling k is irrelevant at 
the SC+SDW to SC critical point. To see this, imagine that the critical point is 
described by the k = point, where the and ^i )2 fields are decoupled. This 
point is actually a fixed point of the renormalization group and is invariant 
under scaling transformations. From (44) we see that the scaling dimensions 
of ^2 is 1, while that of l^a] 2 is (3 — 1/V), where v is the correlation length 
exponent. So from (45) we conclude that the scaling dimension of k is (1/ v— 2). 
The exponent v is determined by a theory like that in (9) for \I/ a and all such 
theories have v > 1/2, the mean-field exponent. Consequently the scaling 
dimension of k is negative, and the decoupled k = fixed point describes the 
SC+SDW to SC quantum phase transition. The critical magnetic fluctuations 
are associated with $ a , and the fermionic, S = 1/2 quasiparticles ^1,2 are 
merely innocent spectators. 

This last result is quite powerful, as it implies that we can apply the results 
of the quantum transition in the insulator, discussed in Section 2, essentially 
unchanged to the magnetic fluctuations in the cuprate high temperature su- 
perconductors. Such a proposal was first made in Ref. [23,24], and leads to 
the prediction that the S = 1 exciton should survive as a stable excitation in 
the SC state near its transition to the SC+SDW state, and that this should 
be evident as a sharp "resonance" in dynamic neutron scattering: such a reso- 
nance appears to have been observed in a number of experiments [25-29] . The 
framework of the SC+SDW to SC transition has also been used to understand 
dynamic spin correlations in the presence of Zn impurities [30,31], and to ex- 




(45) 
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plore the phase diagram in the presence of an applied magnetic field [32,33]: 
the results have been quantitatively compared with neutron scattering [34-38] 
and scanning tunnelling microscopy experiments [39]. 



6 Transition between BCS superconductors 

In this final section, we will consider one of the simplest examples of a non- 
trivial quantum critical point in two spatial dimensions which, surprisingly, 
was only discovered recently [41,42]. The two phases on either side of the 
critical point are described by an entirely conventional BCS theory: they differ 
only in the symmetry of the Cooper pair wavefunction, dependent upon the 
relative co-ordinates of the paired electrons. At T = 0, BCS theory provides 
an accurate description of the low energy properties everywhere except at the 
point which separates the two phases: this is probably the only known example 
of the failure of BCS theory in two (or higher) dimensions in a superconducting 
ground state. At T > 0, this failure broadens into the "quantum critical" 
region which should be easily experimentally detectable: we will not discuss 
these T > issues here, and instead refer the reader to other reviews [1]. 

Our original motivation for examining transitions between BCS supercon- 
ductors was provided by photoemission experiments [40] which observed an 
anomalous broadening of the quasiparticle spectrum. We will discuss this fur- 
ther in Section 6.3, but will motivate the problem now on theoretical grounds. 
The Hubbard model on a square lattice is believed to provide a good start- 
ing point for studying the physics of the cuprate superconductors. There is 
now good evidence that for electron densities close to half-filling, the primary 
instability of the Fermi liquid state is to a BCS superconductor with d x i_ y i 
symmetry [43], whose mean-field theory was presented in Section 5. The op- 
posite limit of a low density of electrons, well away from half-filling, has also 
been examined theoretically [44,45], and the low density expansion permits 
reliable predictions to be made: it was found that the ground state was again 
a BCS superconductor, but now the pairs had a distinct d xy symmetry (in 
the presence of a square lattice, the d x 2_ y 2 and d xy states are not degener- 
ate, as they would be in free space). Here we will answer the simple question: 
how does the ground state evolve from between these limits with qualitatively 
different wavefunctions ? 

Rather than working the full complexities of the Hubbard model, we follow the 
strategy of Section 5 and work with the simplest phenomenological model that 
contains the two phases we are interested in. We want to extend H t j in (38) so 
that a BCS mean-field theory will permit a region with d xy superconductivity. 
Clearly this will be possible if we have a pairing interaction along the (1, ±1) 
directions of the square lattice; so we extend H tJ to allow for a second-neighbor 
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Fig. 11. Values of the pairing amplitudes, —{c^Cji — Cj|Cjf) with i the central site, 
and j is one of its 8 near neighbors. 

exchange along the diagonals, J 2 , as in (23): 

HtJ = E £kci a c k(T + Ji E s i • s j + J 2 E S * " S i- ( 46 ) 

fe (ij) nnn ij 



We will follow the evolution of the ground state of H t j as a function of J2/ Ji. 
We expect that the universal properties of this route will be similar to that 
of the Hubbard model as a function of carrier concentration, in light of the 
studies of the latter model quoted above. 



6.1 BCS Theory 



We begin with a standard BCS analysis of H t j which closely parallels that 
presented in Section 5. The mean-field Hamiltonian is now modified from (40) 
to 

Hbcs = E £ kc\ a c ka - y E A ^( c jt c ]+/U ~ C ji c i+A,t) + h - c - 

k j,n 

~ yE'A^Vj^ - c}^) + h.c, (47) 

where the second summation over v is along the diagonal neighbors x + y 
and —x + y. To obtain d xy pairing along the diagonals, we choose A x+y = 
—A_ x+y = A xy . We summarize our choices for the spatial structure of the 
pairing amplitudes (which determine the Cooper pair wavefunction) in Fig 11. 
The values of A x 2_ y 2 and A xy are to be determined by minimizing the ground 
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Fig. 12. BCS solution of the phenomenological Hamiltonian Htj in (46). Shown are 
the optimum values of the pairing amplitudes |A a ,2_„2| and |A XJ/ | as a function of J 2 
for t\ = 1, t 2 = —0.25, n = —1.25, and J\ fixed at J\ = 0.4. The relative phase of the 
pairing amplitudes was always found to obey (50). The dashed lines denote locations 
of phase transitions between d x 2_ y 2, d x i_ y i + id xy , and d xy superconductors. The 
pairing amplitudes vanish linearly at the two critical points, corresponding to the 
exponent Pbcs = 1 in (53); the slight rounding at the end points apparent in the 
figure is a consequence of the finite step-size used in numerical evaluation of the 
integral in (48). Combining these results with those of [43-45] we propose that 
the more realistic Hubbard model will show the same sequence of phases as those 
above as a function of increasing doping. Dagan and Deutscher [46] were the first to 
propose that d x 2_ y 2 +id xy superconductivity occurs at a doping larger than d x 2_ y 2 
superconductivity on the basis of experiments we will discuss in Section 6.3. 

state energy (generalizing (41)) 

f d 2 k 

Ebcs = Ji\A x 2_ y 2\ + J 2 \A xy \ — J [E k — e k ] (48) 



where the quasiparticle dispersion is now (generalizing (42)) 



Eh 



e\ + I Ji A x 2_ y 2 (cos k x — cos k y ) + 2J 2 A xy sin k x sin k y 



1/2 



(49) 



Notice that the energy depends upon the relative phase of A x 2_ y 2 and A x , 
this phase is therefore an observable property of the ground state. 



xy 



It is a simple matter to numerically carry out the minimization of (49), and 
the results for a typical choice of parameters are shown in Fig 12 as a function 
J2/J1. One of the two amplitudes A 

ground state is always superconducting. The transition from pure d. 
conductivity to pure d xy superconductivity occurs via an intermediate phase 



y 2 or A xy is always non-zero and so the 

y 2 super- 
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in which both order parameters are non-zero. Furthermore, in this regime, their 
relative phase is found to be pinned to ±7r/2 i.e. 

arg(A^) = arg(A ;E 2_j /2 ) ± tt/2 (50) 



The reason for this pinning can be intuitively seen from (49): only for these 
values of the relative phase does the equation Ek = never have a solution. 
In other words, the gapless nodal quasiparticles of the d x 2_ y 2 superconductor 
(and also those of the d xy superconductor) acquire a finite energy gap when 
a secondary pairing with relative phase ±7r/2 develops. By a level repulsion 
picture, we can expect that gapping out the low energy excitations should help 
lower the energy of the ground state. The intermediate phase obeying (50) is 
called a d x 2_ y 2 + id xy superconductor. 

The choice of the sign in (50) leads to an overall two-fold degeneracy in the 
choice of the wavefunction for the d x 2_ y 2 + id xy superconductor. This choice 
is related to the breaking of time-reversal symmetry, and implies that the 
d x 2_ y 2 + id xy phase is characterized by the non-zero expectation value of a Z2 
Ising order parameter; the expectation value of this order vanishes in the two 
phases (the d x 2_ y 2 and d xy superconductors) on either side of the d x 2„ y 2 +id xy 
superconductor. As is conventional, we will represent the Ising order by a 
real scalar field 0. Fluctuations of become critical near both of the phase 
boundaries in Fig 12. As we will explain in Section 6.2 below, the critical 
theory for fluctuations is not the usual 4 field theory which describes the 
ordinary Ising transition in three spacetime dimensions. 

Near the phase boundary from d x 2_ y 2 to d x 2_ y 2 + id xy superconductivity it is 
clear that we can identify 

= i&xy, (51) 



(in the gauge where A x 2_ y 2 is real). We can now expand Ebcs in (48) for 
small (with A x 2_ y 2 finite) and find a series with the structure [47,48] 



£ B cs = £o + s0 2 + t;|0| 3 + ..., (52) 



where s, v are coefficients and the ellipses represent regular higher order terms 
in even powers of 0; s can have either sign, whereas v is always positive. 
Notice the non-analytic |0| 3 term that appears in the BCS theory — this 
arises from an infrared singularity in the integral in (48) over E k at the four 
nodal points of the d x 2_ y 2 superconductor, and is a preliminary indication that 
the transition differs from that in the ordinary Ising model. We will achieve 
a deeper understanding of this non-analyticity when we consider fluctuations 
in Section 6.2. We can optimize by minimizing E BC s in (52) — this shows 
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Fig. 13. Conjectured T = phase diagram for the cuprate superconductors as a 
function of increasing hole concentration 5. The theory for the SC+SDW to SDW 
transition appears in Section 5, while that of the transition describing the onset 
of d x 2_ y 2 + id xy order appears in Section 6.2. The second transition may not be 
on the physical axis for a given cuprate superconductor, but may be nearby in a 
generalized parameter space. 

that ((p) = for s > 0, and (0) ^ for s < 0. So s ~ {JijJ\) c ~ J2/J1 where 
(J 2 /Ji) c is the first critical value in Fig 12. Near this critical point, we find 



where we have allowed for the fact that fluctuation corrections will shift the 
critical point from s = to s = s c . The present BCS theory yields the exponent 
Pbcs — lj this differs from the usual mean- field exponent (3mf — 1/2, and 
this is of course due to the non-analytic |</>| 3 term in (52). 

An essentially identical structure appears at the second critical point in Fig 12 
at the boundary between the d x 2_ y 2+id xy and d xy superconductors. The d x 2_ y 2 
and d xy pairing amplitudes exchange roles, and the Ising order parameter for 
the transition is (f) = iA x 2_ y 2 (replacing (51)). We will therefore not discuss 
this case further. 

We close this subsection by presenting a unified overview of the quantum 
phase transitions in the cuprate superconductors we have introduced so far: 
this is presented in Fig 13. At low, but finite, doping we have the SC+SDW 
phase which was discussed in Section 5. Here the presence of the SDW order 
induces spatial anisotropy between the x and y axis (in real space), and so 
the superconductivity is not purely d-wave, although gapless nodal fermions 
may be present [49,50]. With increasing 5 we have a transition to a SC phase, 
where the superconducting order is d x 2_ y 2, and is described by a vanilla BCS 
theory: this transition was also discussed in Section 5. Finally, the results 
of [44,45] show that the d xy superconductivity at very low electron density: we 
have used a phenomenological model of this crossover in the present section 
to argue that this implies that the next phase with increasing doping should 
be a d x 2_ y 2 + id xy superconductor. The full critical theory of this transition 
appears in the following subsection. In placing the d x 2_ y 2+id xy at large doping, 
we are also following the proposal and experimental results of Dagan and 
Deutscher [46] which we will discuss further in Section 6.3. 



(4>) ~ (Sc - Sf, 



(53) 
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6.2 Quantum field theory 



We will now present a full theory for the quantum transition between the 
d x 2_ y 2 and d x i_ y i + id xy superconductors: this will allow us to compute the 
corrections to the predictions of the BCS theory implied by (52). 



We have laid much of the ground work for the required field theory in Section 5. 
In addition to the order parameter 0, the field theory should also involve the 
low energy nodal fermions of the d x 2_ y 2 superconductor, as described by 
in (44). For the fluctuations, we write down the usual terms permitted near 
a phase transition with Ising symmetry, and similar to those in (9): 



d rdr 



((«9 r 0) 2 + c 2 {d x <P) 2 + c 2 (^0) 2 + S 2 ) + 



u 

24 ( 



(54) 



Note that, unlike (52), we do not have any non-analytic |0| 3 terms in the ac- 
tion: this is because we have not integrated out the low energy nodal fermions, 
and the terms in (54) are viewed as arising from high energy fermions away 
from the nodal points. Finally, we need to couple the and ^2 excitations. 
Their coupling is already contained in the last term in (47): expressing this in 
terms of the fermions using (43) we obtain 

= J d 2 rdr [A0 (^\t v ^ 1 - vIt v V 2 )] , (55) 

where A is a coupling constant. The partition function of the full theory is 
now 

Z did = J V < j)Vm 1 Vm 2 exp -S4- S^) , (56) 



where was in (44). It can now be checked that if we integrate out the ^2 
fermions for a spacetime independent 0, we do indeed obtain a |0| 3 term in 
the effective potential for 0. 

We begin our analysis of Z did in (56) by following the procedure discussed 
below (45). Assume that the transition is described by a fixed point with A = 0: 
then as in Section 5, the theory for the transition would be the ordinary 4 
field theory S^, and the nodal fermions would again be innocent spectators (we 
already know from the presence of the |0| 3 term in (52) that this assumption 
cannot be correct, but we wish to see how it is invalidated in a renormalization 
group analysis). The scaling dimension of at such a fixed point is (1 + r]i)/2 
(where r]j is the anomalous order parameter exponent at the critical point of 
the ordinary three dimensional Ising model), while that of V&i^ is 1 (as before). 
Consequently scaling dimension of A is (1 — rj I )/2 > (using the known, very 
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small value of t]i). This positive scaling dimension implies that A is relevant 
and the A = fixed point is unstable. 



Determining the correct critical behavior now requires a full renormalization 
group analysis of Z did . This has been described in some detail in [51], and we 
will not reproduce the details here. The main result we need for our purposes is 
that couplings A, u, vf/c and i>a/c all reach non-zero fixed point values which 
define a critical point in a new universality class. These fixed point values, 
and the corresponding critical exponents, can be determined in expansions in 
either (3 — d) [41,42] (where d is the spatial dimensionality) or 1/N [52] (where 
N is the number of fermion species). Indeed the fixed point has the structure 
of the so-called Higgs- Yukawa model which has been studied extensively in the 
particle physics literature [53] in a different physical context: quantum Monte 
Carlo simulation of this model also exist [54], and provide probably the most 
accurate estimate of the exponents. 

The implications of the existence of this finite-coupling critical point are sim- 
ilar to those for (9) in (13). The fermion correlation function G\ = 

obeys 



Gl(fc,W, "(^ + ^-^ )/2 ( j 



at low frequencies for s > s c . Away from the critical point in the d x 2_ y 2 
superconductor with s > s c , (57) holds with rjf — 0, and this is the BCS 
result, with sharp quasi-particle poles in the Green's function. At the critical 
point s = s c (57) holds with the fixed point values for the velocities (which 
satisfy vp — v& — c) and with the anomalous dimension i]f ^ — the (3 — d) 
expansion [41] estimate is rjf ~ (3— d)/14, and the 1/N expansion estimate [52] 
is rjf pa l/(37r 2 iV), with N = 2. This is clearly non-BCS behavior, and the 
fermionic quasiparticle pole in the spectral function has been replaced by a 
branch-cut representing the continuum of critical excitations. The corrections 
to BCS extend also to correlations of the Ising order 0: its expectation value 
vanishes as (53) with the Monte Carlo estimate (3 ~ 0.877 [54]. The critical 
point correlators of also obey (13) with the exponent rj ps 0.754 [54], which is 
clearly different from the very small value of the exponent rjj at the unstable 
A = fixed point. The value of (3 is related to rj by the usual scaling law 
/3 — (1 + 77)^/2, with v ~ 1.00 the correlation length exponent (which also 
differs from the exponent vj of the Ising model). 
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6.3 Connection with experiments 

We noted at the beginning of this section that our original motivation for 
studying the d x i_ y i to d x 2_ y 2 +id xy transition was provided by photoemission 
experiments of Valla et al. [40]. In particular these experiments found that 
the typical quasiparticle energy width was of order ksT. As we have reviewed 
elsewhere [1], this is a characteristic property of the T > quantum critical 
region of a T = critical point which obeys strong hyperscaling properties. In 
particular, it is essential that the critical point be described by an interacting 
field theory, with the interaction strengths determined by a finite-coupling 
fixed point of the renormalization group. This condition is clearly satisfied by 
Zdid'- the breakdown of BCS theory at s = s c , which leads to (57) with r\f ^ 
at T = 0, is directly responsible for quasiparticle lifetimes of order h/ksT over 
a finite range of g values in the T > quantum critical region. In Ref. [41] 
it was also shown that the fixed point of Z did was also essentially unique in 
satisfying such requirements: the only other zero momentum order parameter 
which led to suitable critical point was that between a d x 2_ y 2 and a d x 2_ y 2 +is 
superconductor. 

A more direct probe of the d x 2_ y 2+id xy order has been provided by very inter- 
esting tunneling measurements of Dagan and Deutscher [46]. They examined 
the splitting of a zero-bias conductance peak as a function of doping S and an 
applied magnetic field, H, perpendicular to the two-dimensional electron gas. 
They found that this splitting closely tracked the value (0) as determined by 
minimizing the energy in (52) extended with a linear coupling to the magnetic 
field 

Ebcs = E + s<j) 2 + v\<j)\ 3 - H<t>. (58) 



The linear coupling is also just that expected from the fact that the d x 2_ y 2 + 
id xy superconductor has a spontaneous orbital magnetic moment [55,48]. The 
experiments of Ref. [46] are able to scan across the critical point at s = 0, and 
observe behavior consistent with (58); we regard this as strong evidence for 
the phase diagram of Fig 13 with a quantum critical point between d x 2_ y 2 and 
d x 2_ y 2+id xy superconductors. It would be interesting if future experiments are 
able to test for the corrections to (58) discussed in Section 6.2. We also note 
here the recent measurements of Yeh et al. [56] presenting evidence for related 
(but somewhat different) changes in the pairing symmetry in the overdoped 
regime. 
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7 Conclusions 



We have studied a number of phase transitions of two dimensional correlated 
electron systems and a schematic phase diagram summarizing their relevance 
to the high temperature superconductors was in Fig 13. At zero doping, S = 0, 
we start from an insulating antiferromagnet with long range Neel order. After 
a localized carrier regime at very low doping, the ground state at small S 
appears to be superconducting with co-existing collinear spin-density wave 
order (the SC+SDW phase). There is a long in-plane spin correlation length in 
this phase (the finiteness of the correlation length is, we believe, due to disorder 
effects, and so at long enough scales the spin order is spin-glass like). At special 
commensurate dopings (like 5 — 1/8) the inter-plane coupling is strong enough 
to induce three-dimensional long-range order at finite temperatures. 

At larger 5 there is a transition to a pure d x 2_ y 2 superconductor, and this 
transition was described in Section 5. Charge order and its fluctuations also 
play a central role in the magnetic ordering quantum transition: we discussed 
these effects in the insulator in Section 4. They also have an important effect 
near the corresponding SC+SDW to SC transition in the superconductor, but 
the reader is referred to [21] for a review of these. 

With increasing 5, we argued on the basis of theoretical [44,45] and experimen- 
tal studies [46] that the cuprate superconductors are at least in the vicinity 
of a separate transition from the d x 2_ y 2 superconductor to a d x 2_ y 2 + id xy su- 
perconductor. A quantum field theory for this transition, associated with the 
breakdown of BCS theory at the critical point, was presented in Section 6. 
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